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to Maxwell's equations. The results are an ideal starting point for calculations of 
corrections to the paraxial approximation. 
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1. Introduction 

Solutions to tlie wave equation in spheroidal coordinates have application to a wide range 
of problems in physics [1]. In particular, in the short wavelength limit, consideration 
of the wave equation in oblate spheroidal coordinates leads directly to the well 
known Gauss-Laguerre solutions to the paraxial approximation of the wave equation 
in optics [2]. The paraxial approximation is an extremely versatile tool for studying 
beams of coherent radiation (e.g. laser beams), and successfully describes phenomena 
such as propagation of beams through lens systems, focus size, wave front curvature 
and phase shifts in the focus of beams, as well as eigenmodes and eigenfrequencies of 
spherical mirror resonators [3]. Despite its great success, the paraxial approximation is 
only an approximation. Consideration of corrections beyond the paraxial approximation 
is relevant not only in order to establish a bound on its validity but is in fact necessary 
to describe strongly focused beams [1] and to explain the fine structure of the frequency 
spectrum of high finesse resonators [5] . 

Corrections to the paraxial approximation can be calculated by reconsidering the 
term which is neglected when deriving the paraxial approximation from the wave 
equation [6]. As an alternative, corrections can be effectively considered using exact 
solutions to the wave equation in spheroidal coordinates. This has the advantage 
of removing ambiguities in the definition of higher-order terms, leads to simpler 
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expressions, and is in general a more natural framework for exact consideration of 
beamlike solutions to the wave equation. 

The wave equation is separable in spheroidal coordinates, allowing solutions to be 
written as the product of a so called radial function depending only on the ^ coordinate, 
a so called angle function depending only on the rj coordinate and the function e™**^ 
depending only on the azimuthal (j) coordinate [1]. For short wavelengths, solutions for 
the radial and angle functions in the form of asymptotic expansions have been known 
for a long time [71 El E]- The close relationship between spheroidal coordinates and 
Laguerre- Gaussian beams suggests that the asymptotic solutions are in some way related 
to the Gauss-Laguerre solutions to the paraxial wave equation. While this relationship 
is directly obvious for the asymptotic solutions for the angle functions, the opposite is 
true for the previously known solutions for the radial functions. In fact, the previously 
known asymptotic solutions for the radial functions consist of multiple sums with a huge 
number of terms. 

Using solutions to the wave equation for exact calculations in optics requires taking 
the transverse vector character of the electromagnetic fields into account. Full vector 
solutions to Maxwell's equations were already considered by Ref. ^U\. These vector 
solutions are derived by applying vector operators to the scalar spheroidal solutions and 
are expressed as derivatives in the spheroidal coordinates of the scalar solutions. As 
for the asymptotic solutions for the radial functions, the resulting expressions are quite 
complicated. 

In this paper, we derive a new expression for the asymptotic expansion of the 
radial functions which to lowest order is directly reminiscent of solutions of the 
paraxial approximation. The new expression is significantly simpler, making its analytic 
application practicable. In addition, we derive expressions for the Cartesian derivatives 
of the spheroidal wave functions which allow us to define a significantly simpler set of 
vector spheroidal wave functions. The component of the field transverse to the direction 
of propagation is equal to a single scalar spheroidal wave function. This creates a 
direct correspondence between our vector solutions and Gauss-Laguerre solutions to the 
paraxial approximation. 

After a brief discussion on notation used in this paper in section [TTTl we introduce 
oblate spheroidal coordinates in section [21 The wave equation in spheroidal coordinates 
as well as separation of coordinates is discussed in section 12. 1[ Section [3] is a review 
of the asymptotic expansion of the angle functions. The derivation is equivalent to the 
one given e.g. in Ref. [1]. The radial functions are considered in section HI Starting 
with the previously known asymptotic expansion allows us to derive a new expression 
for the two lowest-order terms. Taking the lowest-order term as an ansatz for the radial 
function leads to a differential equation which allows higher-order terms to be calculated 
iteratively. In section [5] we define scalar spheroidal wave functions and demonstrate their 
equivalence in the short wavelength limit to the Gauss-Laguerre solutions to the paraxial 
wave equation. This is followed by a brief discussion concerning the symmetry properties 
of the scalar spheroidal wave functions upon reflection through a plane containing the 
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symmetry axis of the spheroidal coordinate system. We proceed in section O by deriving 
expressions for the Cartesian derivatives of the scalar spheroidal wave functions. These 
expressions allow us to define a new set of vector spheroidal wave functions in section [7] 
which are transverse vector solutions to the wave equation, section [7] concludes with 
expressions for the curl of the vector spheroidal wave functions. 

1.1. Notation 

The notation used in this paper is generally consistent with the notation in Ref. [1], 
with the exceptions noted here. Previously, the symbol c has been used to quantify the 
scaling of the spheroidal coordinate system relative to the wavelength. Since it seems 
problematic to use the same symbol as for the speed of light in a theory with strong 
application to optics, we use c instead. 

Additionally, Ref. [1] generally uses m and n to label the spheroidal functions 
and associated parameters, using the labels m and u with (m, n) = (m, 2z/ + m) only 
when Laguerre functions are directly involved. In the short wavelength limit, both sets 
of labels can be used. The advantage of the m, u labeling is that it is identical to 
the labeling of the Laguerre polynomials, which are essential for the expansion of the 
angle functions in the short wavelength limit. As a result, we use the labels m and u 
exclusively. 

Outside the short wavelength limit, not considered in this paper, use of the labels 
m and n is necessary. This is due to the fact that when the short wavelength limit does 
not apply, the labels m, v do not uniquely label the set of solutions to the wave equation 
whereas the labels m, n do. Each pair of values m, v corresponds to two solutions for 
the angle function and two solutions for the radial function with opposite parity about 
the plane ^ = 0. In the short wavelength limit, the two distinct angle functions converge 
to the same function asymptotically, and the two radial functions can be considered as 
the real and imaginary part of a single complex valued function, so that m, v effectively 
labels all solutions uniquely. 

In contrast to Ref. [1], we attempt to consistently place the labels m and v as 
subscripts, generally placing other labels as superscripts. The only exception is the 
Laguerre functions, where we use the standard notation lI™"*. Since we are only 
interested in oblate spheroidal coordinates and real arguments r/ and we write the 
angle and radial functions respectively as Smvij]) and RmuiO instead of as Smn{—ic,ri) 
and Rmn{—ic,iC,) as in Ref. [1], thereby ignoring the relationship between solutions in 
oblate spheroidal coordinates and in prolate spheroidal coordinates. Finally, the indices 
for the sums over Laguerre functions as e.g. in Eq. f|T2l) range from to cxo rather than 
from —u to oo as in Ref. [1]. This allows us to treat all solutions for the angle functions 
on an equal footing, without choosing a specific one in advance by fixing u ahead of 
time. 



Solutions to Maxwell's Equations using Spheroidal Coordinates 



4 



2. Spheroidal coordinates 

The two-dimensional elliptic coor- 
dinate system is defined from the 
set of all ellipses and all hyperbo- 
las with a common set of two fo- 
cal points [1]. We denote the sep- 
aration of the two focal points by 
d. Oblate spheroidal coordinates 
are derived from elliptic coordinates 
by rotating the elliptical coordinate 
system about the perpendicular bi- 
sector of the focal points [1]. The 
focal points thereby sweep out a cir- 
cle of radius (i/2. Oblate spheroidal 
coordinates are commonly mapped 
to Cartesian (x, y, z) coordinates by 
placing this circle in the x — y plane 
with center at the origin. The co- 
ordinates are often labeled ?7, ^ and 
with the transformation to Carte- 
sian coordinates given by [T] 

d 




Figure 1. Oblate Spheroidal Coordinates. 
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Oblate spheroidal coordinates reduced to the x 
The coordinates r and 6 defined by 

r = \/x^ 



and 



cos 6' 



z plane are shown in Fig. ([T]) 



V 



(2) 



are useful to describe aspects of the coordinate system. Note that near the 2;-axis, ^ 
quantifies distance along the z-axis and 6 quantifies distance from the 2;-axis. To lowest 



order in 9, z = and r = ^^/TT~^0. 



2.1. Wave equation in spheroidal coordinates 

The wave equation, 

(V^ + p)^ = 0, 



(3) 
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written in oblate spheroidal coordinates, is given by pQ 

(4) 



_ kd 
with c = — . 

2 

This equation is separable. Therefore ip can be written as a product of three functions 
depending only on r], ^ and 0, respectively. Due to cylindrical symmetry, the function 
depending on (p can be written as e^"^"^ for integer values of m, 

= Rrnu{OS,nu{vy'^^- (5) 

The label v is used to denote the various solutions of Eq. (jlj). 

The so called radial function RmuiO angle function Smuijf) satisfy the two 
differential equations 

+ ^')| + + <^'^" - ^-(0 = (6) 

and 

as can be seen from comparison with Eq. (jl]). Here, \mv is the separation constant. 



3. Angle functions for large c 



In the limit c — >• cx), a solution to Eq. ([7]) expressed as a series in | can be found as 
follows. The ansatz 



Smiri) = i^-V 
allows Eq. ([7]) to be rewritten as 
d' . 



fe-5(i-'?)s 



with 



X = 2c(l — 77), 



(8) 



X 



dx'^ 



(m 



dx 



m 



1 c^ + A 



4c 



4c 



+ 



X 



— 2x{m + 1) d x{m + 1) — (m^ + m) 



4c 



dx 



+ 



4c 



(9) 



We have dropped the index u in order to refer to a general solution with arbitrary Xm{u). 
Note that by introducing the factor (l—ri'^)~ in Eq. (JHl), we break the symmetry between 
positive and negative m. This symmetry nonetheless persists in the original Eq. ([7]) so 
that the solutions must in the end be symmetric under a transformation from m to —m. 
We address this issue in more detail in section [51 

For c — s> 00, the terms on the second line of Eq. inside the bracket vanish. 
What is left is simply the associated Laguerre differential equation. A possible set of 
approximate solutions to Eq. (Q is therefore the set of associated Laguerre polynomials 



x) with fixed m. 



1, 



(10) 
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and Xm given by 

A™ = -c' + [2(m + l)+4z/]c + C(l). (11) 

The form of Eq. ( |T0|) suggests expanding the exact solution Sm{x) in terms of associated 
Laguerre polynomials as 

oo 

.„(x) = 5^A^L^)(x). (12) 

s=0 

The Af^ are expansion coefficients. Inserting this expression into Eq. iQ and using 
appropriate recursion relations for Laguerre polynomials allows the left hand side of 
Eq. ([9]) to be rewritten as a linear combination of Laguerre polynomials with fixed m. 
Note that a number of useful recursion relations for Laguerre Polynomials are listed in 



Appendix A[ Since the Laguerre polynomials with fixed m are linearly independent, 
the coefficient of each term in this sum must be zero, leading to the following relations 
between the A^, 

s{s + m)yl-i + [(4. + 2(m + l))c - {2s' + (2s + l)(m + 1)) - \'JA:^+ 

{s + l)is + m + l)A:^' = 0, (13) 

with \'^ = \rn + c'- 

These relations can be written in matrix form as 

(cMo + Ml - A;,I)A = 0. (14) 

Here, Mq is a diagonal matrix. Mi is a tridiagonal matrix, I is the identity matrix and 
A is the vector of coefficients according to 

{Mo)s,t=[^s + 2im + l)]Ss,t, 
{Mi)s,t = s{s + m)6s^i,t - (2s2 + (2s + l)(m + l))6s,t + {s + l){s + m + l)Ss+i,t, (15) 
(A). = A^. 

We use the Kronecker delta function, Ss,t- Eq. ( fT4l) transforms the problem of 
determining the coefficients in Eq. (fT2l) into an eigenvalue problem. A^ is an eigenvalue 
of cMq + Ml and A is the corresponding eigenvector. 

In the limit c oo, Mi can be neglected compared to cMq. Since Mq is diagonal, 
finding an approximate solution to Eq. f|T^ is therefore trivial, this solution being given 
by Eqs. ( fTOj) and ( fTTi) . Improved approximate solutions can be found via perturbation 
theory [11]. According to convention we reintroduce the label u to denote the solution 
whose eigenvector is given to lowest order in | by A^^ = dg^u- We expand A^^, and A^^^ 
in powers of 4 as 



k=-\ 
oo 



k --k 



(16) 



s,k z^—k 
mu 



k=0 
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Inserting these expressions into Eq. (I14p and matching terms of equal order in ^ leads 
to recursive expressions for and A^^, 

/?L = X^(MO.,.A:;t k > 0, 

AS,k+l _ l^t\^^^^)s,t^mv Z^k'=0 l-^mu^mu U ~^ n a -I- i, 

"^^ " (Mo)(.,) - (Mo)(.,) k>^.s + v. 

Since M\ is tridiagonal, is zero for |s — z/| > k. As a result, each sum over s or t 
contains only a finite number of nonzero terms. 

Eqs. (ITTl) lead to the expressions for fi^^ and for /c < 4 given explicitly in 
Ref. [T]. Higher-order values can easily be obtained numerically from Eqs. fll7l) . Note 
that the expression for fH^"^ in Ref. [1] contains an incorrect " — " sign. 

4. Radial functions for large c 

Following a similar procedure as the one in section [31 it is also possible to find an 
asymptotic expansion for the radial function, given by [1] 

RmAi) = ^"^t^Mo"^"^'' ^-^"^'^^^ne^ + ir'f:A^m.Ui-^[2-c{l + (18) 

mu g=0 

Note that Ref. [T] defines a total of four different radial functions Rmn{—ic^i^) with 
j = 1,...,4. These are respectively the real part of Eq. f|T8|) . the imaginary part of 
Eq. ffTSl) . Eq. fITS]) itself and the complex conjugate of Eq. (ITS]) . Without loss of 
generality, we only consider Eq. (ITS!) . The A^^^ in Eq. (fTS!) are the expansion coefficients 
for the angle function from Eq. (IT^ which can be calculated using Eqs. fll6p and (I17p . 
The f/i"^^(2) are second solutions of the Laguerre differential equation, defined by [12] 



^-imn J {-m) ( \ _ V{u+m+l) rri j (-m) / ^\ 

^M^^)^^! "^-^"-^'^ (19) 

sin(m7r) 

For integer values of m, the denominator and numerator in Eq. (fT9l) both equal zero, 
in which case u!)^\z) is defined as a limit in m. For large \z\, an asymptotic expansion 
for U^\z) exists, given by pi fT2] 

^ e'(-+V^)-e- - + ^)! + ^ + ^)! 

r=0 

Unlike the asymptotic expansion for the angle function, which is quite transparent 
since it contains a single Laguerre polynomial to lowest order and additional Laguerre 
polynomials as higher-order corrections, the asymptotic expansion for the radial function 
is rather opaque since the terms s = through s = z/ all contain terms of equal and 
lowest order in | as can be seen from Eqs. (ITSl) and (l20ll . We attempt to rewrite Eq. (fTSll 
as a sum over terms of equal order in 4. To this end we insert the asymptotic expansions 
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for u!)^\z) and A^^^, Eqs. ( !20!) and ( fT6|) . into Eq. ( fTSj) . This leads to sums over r, s 
and fc. We define p = r + s + k — u and observe that the total order in | of a term in 
the triple sum is p + 1. Replacing the sum over by a sum over p, one obtains 



m!r!s! A^, (1 + ^0 



m+r+s+l 



The symbol [.J denotes rounding down to the nearest integer. 

As we will demonstrate, the sums over r and s can be performed resulting in 
relatively compact analytic expressions for all values of p. We start by considering the 
two lowest-order cases, p = and p = 1. To proceed, we need analytic expressions for 
A'Jl foT s + k<u + l. Using Eqs. it can be shown by induction that 

[v — sy.sXvm + s ! 



and 



,,,^-,+1 _ ^ 2(u~s+V) (3z/ + 2m + g + V)v\{v + m)! 

(z/-s-l)!s!(m + s)! ' ^ ^ 



For p = 0, one obtains 



V 



m+s+l 



c-(^+i) (z/ + m)! ^.gg (1 - i^)'^W2 



(24) 



-e 



Al^^ 22^m! (1+^0''+"/^+^' 
The first line is Eq. fl2Tl) with fixed p = 0. To obtain the second line, define z = \ ^ i^ 
and observe that the first line is essentially the binomial expansion of (1 — |)''. 
The case p = 1 is already slightly more complicated. 



1 V 



c-(-+2) (y + 



AO,, 22'^m! 



(z/ + l)(z/ + m + 1) 3z^ + 2mz/ + z/ z/(z/ + m 



^1 _ ,^)-^W2 

2(1 + zO ' 4 2(1-^0 J" (1 + z^)^W2+i ■ 

The first line is again Eq. (I2T]) but with fixed p = 1. The second line is obtained in 
a similar manner as before, by observing that for r = the first line is essentially the 
binomial expansion of 2:^''+™'+^^ 2r~(^'^+2m) ^ i^^^ r = 1 the first line is 

essentially the binomial expansion of ^ z''^^^'^^ (l — l)'^ . 

Using Eqs. ([23]) and fl2S]l . we obtain a preliminary result for the summation of 
terms of equal order in | for the asymptotic expansion of RmuiO through first order 
in |. In writing the preliminary result, we must confront the issue that the higher- 
order terms are to some degree arbitrary. RmuiO defined by its differential equation 
only up to some constant factor. In particular, we can multiply RmuiO by a constant 
which depends on c. The asymptotic expansion in | of such a product contains different 
higher-order terms. Specifically, a constant a times the zeroth order term can be added 
to the pth order term as long as a times the nth order term is simultaneously added 
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to the (p + n)th order term, resulting in a different asymptotic expansion. This issue is 
raised by the presence of A^^^ as part of the normahzation constant in Eq. f|T8|) . which 
depends on 4 through first order according to 



^un~2u i^ + m)\ ^ --u~in~2u (3^^ + 2m + l)z/(z/ + m)! ^ qi^--u~2-j 



ml 



4m! 



(26) 



Replacing A'^^ according to Eq. (126!) leads to 



c 



(z/ + l)(z/ + m + 1) z/(z/ + m) 



2(1 + ,() 2(1 - 

This is a new and significantly simpler expression for the asymptotic expansion of the 
radial function through first order in |. 

Continuing the previous procedure for larger values of p becomes increasingly 
cumbersome, not least because analytic expressions for for larger values oi s + k — u 
become increasingly difficult to find. The similar form of the results for p = and p = 1 
suggests an alternative approach. We factor out the ^ dependent terms in front of the 
bracket in Eq. ([2 



(28) 



Inserting this expression into the differential equation for RmuiO leads to a differential 
equation for r^uiOy 

d^ . ,^ ,^ ^ d 



1 + r)-77T7 + t (2c(l + - 2(2z/ + m + 1))— + 



2(z/ + l)(z/ + m + 1) 2u{iy + m) 



1 + i^ 1-i^ 
with \'^^ = \r,u + - [2(m + 1) + 4z/]c. 
We can write r^uiO an asymptotic expansion in 4, 



d^ 



(29) 



oo 



p=0 



(30) 



Eq. fl29p can be transformed into a recursive expression for the r^l(0 



2(1 +e 



'de 



d 



2z(2z/ + m + 1)— + 
at, 



2(i^ + l)(i^ + m + 1) 2z/(z/ + m) 



/ J t-'mv mv \S/ 
fe=0 



(31) 



di. 



The are the expansion coefficients of \mv according to Eq. (ITBi) . Through p = 2 
one obtains 

{u + + m + 1) v{v + m) 



^(0) 



m v VS y 



2(1 + ^0 



2(1 - ^0 ' 



(2) (c\-^ f [v + + 2){v + m + l){u + m + 2) v{v - l)(z/ + m)(z/ + m - 1) 
A4) - o I (l + z^)2 + 



{u + l){u + m + + mi/ — 2u 
l+ii 



(1 - 

m — 1) z/(z/ + m) (z/^ + mu + + 2m + 2) ~ 



(32) 
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Note that the rml(0 only defined up to a constant of integration. This refiects 
the discussion before Eq. (!26|) . Changing the constant of integration for rml{^) is 
equivalent to adding rml{C,) and changing the higher-order functions accordingly as 
discussed above. We stick to the convention resulting in the presence of in the 
normalization coefficient of Eq. ( ITSl) that \im\^\^oofml{C,) = for p> 0. 

Eq. (13T]) provides a systematic method for determining analytic expressions for the 
summation over r and s in Eq. fl2T|) as claimed directly thereafter. The RmuiO defined 
by Eq. (ITSl) and the Rmv{C) defined by Eq. ( l28l) are solutions to the same differential 
equation with the same boundary conditions and are therefore equal, up to a constant 
factor. The asymptotic expansions of both functions are therefore equivalent so that 
the summations over r and s must be obtainable from the rml{0- 

5. Scalar spheroidal wave functions 

Following Eq. ([5]), we define the scalar spheroidal wave functions ijjmu according to 
^U<P, V, = e-f'x{c-/'il - vT^'e-'^-'-'hUv)) x ^ ^^y^^/,^, e^'^ rUO) ■ (33) 

Note the presence of the factor c™/^ which is included because 1 — ~ ^( 1) (^^^ below) 
and therefore c"'/^{l - r/2)™/2 ~ C(l). 

We briefly demonstrate the equality to lowest order in | of Eq. fl33l) to the Gauss- 
Laguerre solutions to the paraxial wave equation. According to Eq. (j2]) we have 
1 — ?7 ~ 6'^/2. Due to the presence of the term e~'^^^~^\ iprnu takes on finite values 
only for 9 = 0{c-^^^). We can therefore make the following substitutions in Eq. fl33|) . 

sUv) = 4"^H2c(i -v)) + o 0) = L^r\ce') + o 

so that 

x[[cU) e V ^(l + 2z^/d)''W2+i^ J U 



2kdr^Y" f 2kdr^ \ (1 - 2z^/rf)-+-/2 \ ^^fi 



d^ + 4zy \d'^ + 4zy (1 + 2i2/rf)'^+"V2+i 



This is the equation for a Gauss-Laguerre beam with confocal parameter d [3]. 

While the differential equations for the spheroidal functions are independent of the 
sign of m, the approach used to find the asymptotic series for the spheroidal functions 
introduces a dependence on the sign of m. As a result, it is by no means obvious that 



(34) 
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there is some symmetry between positive and negative m in Eq. fl33l) . Demonstrating 
this symmetry is in fact nontriviaL The fact that the Laguerre polynomials transform 
from positive to negative m according to [13] 

Li+mi^) = 7 ~ (35) 

suggests that the spheroidal wave functions remain invariant under the transformation 
(m, z/) — i> (— m, 1/ + m), except possibly for a constant factor. While Eqs. fl28H32l) 
involved in the new asymptotic expansion for the radial functions as well as Eqs. ( |T3lfT7j) 
involving the ^^'(^^ and are all seen to remain invariant under this transformation 
(i.e. (m, z/, s) {—fn, u + m,s + m)), complications arise for the angle functions due to 
the factor the asymptotic expansion of Sm^ and due to the dependence of 

Eq. (!35l) on u and m. These two effects in fact cancel as we checked through first order 
in |, leading to 

V, = i-^rj^;^^ (l + ^^^^ + + (^) ) V, 0- (36) 

Surprisingly the transformation from positive to negative m involves a factor which 
depends on c. 

6. From scalar to vector wave functions 

The scalar spheroidal wave functions defined by Eq. fl55]) are solutions to the wave 
equation. Nonetheless, they are not suitable for exact calculations in optics, since 
solutions to Maxwell's equations are vector functions. Let E be the electric field of 
a solution to Maxwell's equations. In order to specify E, it is useful to choose a vector 
basis and to specify the components of E in this basis. Since our wave functions are 
defined in spheroidal coordinates, it seems reasonable to express E in terms of the 
spheroidal unit vectors, E = E^e^ + -E^e^ + E^e^. Doing so has the major disadvantage 
that the individual components E^, Ej^ and E^ are not solutions to the wave equation 
and therefore have no obvious relationship to the scalar spheroidal wave functions. We 
avoid this problem by using Cartesian unit vectors and write E = i?^.x + Eyj + EzZ. 

A sufficient condition for E to be the electric field of a solution to Maxwell's 
equations in free space is that each of the components E^, Ey and Ez satisfies the 
scalar wave equation and that V ■ E = [H]. As a result, we can construct solutions for 
E from our scalar spheroidal wave functions by writing each of E^., Ey and E^ as a linear 
combination of the iljmv such that V ■ E = is satisfied. Since V ■ E = + + 
satisifying V ■ E = requires knowledge of ^^"^ , ^g™" and ^g"" . Since Cartesian 
derivatives commute with the wave equation, ^g"" , '^g^" and ^g^" are again solutions 
to the wave equation, and can therefore be written as linear combinations of the il)mv 
We try to determine the coefficients of these linear combinations. 
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We begin by deriving expressions for ^ and ^ in spheroidal coordinates. The 
scale factors in spheroidal coordinates, defined by h^^ 



12 



9x 



dui 



h. 



d v^ + e 



hr 



1+e ' 

■oidal 

Cartesian components, are 



1—7] 



2 ' 



with Ui G {^, 7], 0}, are [T] 

(37) 



The spheroidal unit vectors, given by e^^ = t~-§^ and expressed in terms of their 



-irr^ cos 0,^1 



1+e 



r/2 + e 



2 / ' 



1+e 



2 COS(/), -7]^ 



1+e . , , 



1 



(38) 



^2+^2 



§0 = (— sin 0, cos 0, 0) . 
The gradient of a function expressed in spheroidal coordinates is 



1 + 

9a; ' dy 



1 — 7^2 ^ 1 



(9-0 



Obtaining expressions for ^ and ^ is simply a matter of combining Eqs. (I38l) and 



and 



7] d 



± 



9 



(9 _ 2 /^r/(l + e^) d ^ ^(1 -r/2) 9 



(41) 

dz (i \ 772 _|_ ^2 j^2 _|_ ^2 J 

From here on, one can obtain expressions for the Cartesian derivatives of the ipmu 
essentially by brute force. Eqs. (140|) and (14T!) are applied to Eq. (133|) . The result is 
compared to successive orders in 4 to Eq. ( |33l) . One obtains. 



9x dy 



m+l,!/— 1 



l/z/ + m u u + 7n + 1 Z/ + 1 

^ Wm+l,u-2 ~ ^ Wm+l,u-l ^ Wm+l,v H ^ Wm+l,u+l 



1 / (i/ + m)(i/ + m — 1) 



c2 



16 



m+1,1/— 3 



4 r,,.^.,.. ^ 

{u + m)(3z/ + m — 1) 



'0m+l,i/-2 + 



(2z/^ + 4mz/2 + 7z/2 + 2m2i^ + 4mz/ + 3u - 

16 

(2z/^ + 2mi/^ — z/^ — 6mz/ — 5z/ — 2m^ — 5m — 2) 

16 



(42) 



m+l,u—l' 



m+l, I' " 



[z/ + l)(3z/ + 2m + 4) (z/ + l)(z/ + 2) 

7: Wm+l,u+l H 77 '0m+l,i/+2, 
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dx dy ^ ^' 



d 



{U + m)ilm-l,v + (l^ + l)i)m-l,u+l 



1 f {p + m){i' + m — 1] 



[v + m) 



m—lM 



(z/+ l)(z/ + 2) 



1 /(z/ + m)(i/ + m — + m — 2) 
16 



m—l,u—2 



viy + m)(z/ + m — 1) (z/ + m)(2z/'^ + 2mz/^ + z/^ — 2mz/ + 3z/ — m^) 

Wm-l,v-l H 77 Wm-l,v- 



4 16 

{u + l)(2z/^ + 4mi/^ + 5z/^ + 2m^z/ + Amu + 7z/ + 3m + 4) 

16 



m—l,u+l ' 



(z/ + l)(z/ + 2)(z/ + m + 1) (z/ + l)(z/ + 2)(z/ + 3) , 

^ '?/'m-l,i/+2 H 77 ^m-l,i/+3 



and 



d_ _2i 
oz d 



u + m 



16 

2z/ + m + 1 



Z/+ 1 



2 " 2 

1 / (z/ + m)(z/ + m — 1) z^(z/ + m) 6z/^ + 6mz/ + 6z/ + + 3m + 2 

-^m,i.-2 ^ 1pm,u-l 5 ^71 



4 

[z/ + l)(z/ + m + l) (z/ + l)(i/ + 2) 



4 



oil 



These expressions, Eqs. ( I42ll44l) . are extremely powerful. They allow derivatives of the 
scalar spheroidal wave functions to be evaluated with complete disregard to the actual 
structure of the ipmu- As an example, using Eqs. ( I42ll44l) . one can easily check that ipmu 
indeed satisfies the wave equation. 



7. Vector spheroidal wave functions 

Defining a set of vector spheroidal wave functions in terms of their Cartesian components 
is now relatively simple. In analogy to the Gauss-Laguerre solutions to the paraxial 
approximation being the transverse component of the electric field, we set and Ey 
as being proportional to a single iprnu and determine E^ using Eqs. fl42ti44p such that 
V ■ E = is satisfied. Due to their e*'"'^ 0-dependence, the ipmu are eigenfunctions of 
the angular momentum operator about the 2;-axis. We therefore attempt to define the 
vector functions so that they are also eigenfunctions of the angular momentum operator 
about the z-axis. This is accomplished by choosing the transverse vector character of 
the vector wave functions as x + iy and x — iy, which corresponds to cr^ and a~ light, 
respectively. We obtain 

zz 

1 /z. + J-1 , 7U + 4J-2 , 



l/jj,u-l + ^Ju + 



c\ 4 ■ 4 

7z^ + 3 J + 2 , z/ + 1 \ ^^fl 

i'ju + —^i^j,u+i + o 
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E 

1 

c 



Ja V 



(x - iy)ijj+i^^ + 
J)(z/ + J + 1 



4 



J- 
-J 



' l)(5z/ 



J 



3), 



.J,u-l 



(z/+ l)(5z/ + 4J + 7) 



4 

(z/ + l)(z/ + 2) 



o 



(46) 



We use J as the total angular momentum, equal to the sum of the orbital angular 
momentum and the spin angular momentum. 

So far we have ignored the fact that our solutions to the wave equation are 
traveling waves which can travel in both the +^ and direction. All solutions so 
far have consistently represented waves traveling in the +^ direction when an e"*"^* time 
dependence is assumed. This is acknowledged in Eqs. (H5ll and ( H6l) by the + superscript 



label in Ej^^^. We can also define vector spheroidal wave functions Ej^^^ traveling in 
the — ^ direction. These are given by the same expressions, Eqs. ( l45l) and ( H6l) . except 
that z is replaced by — z and the scalar spheroidal wave functions are evaluated at — ^ 
instead of ^. 

For many problems in physics, knowledge of both the electric as well as the magnetic 
field of an electromagnetic beam is needed [H]. Obtaining the magnetic field from the 
electric field requires taking the curl of the electric field and vice versa. The curl of the 
vector spheroidal wave functions can be obtained rather elegantly using Eqs. fl42ti44p to 
evaluate the necessary derivatives. We present the result here as an example for the 
usefulness of these equations. 

For an electric field divided into components according to 



E = (x + iy)E, 

the curl is given by 
'ldE+ 



(x - iy)E^ + zE, 



(47) 



V X E 



I 

I dz ^2 
d 



d_ 

dx 
d 



.d_ 

dy 



E. 



dx dy 



— + i—]E.+-—-i—]E 



X 

d 



d 



.dE_ 

dz 



1 

Yi 



d_ 

dx 



.d_ 
dy 



E., 



dx dy 



Inserting the necessary expressions from Eqs. (1421146 P and writing the result in terms of 



the Ej^^ we obtain 



V X Et , = 

(z/ + J - l)(z/ + J-2) 



E^ 



J -l){2u + J -I) 



E 



Ja+u 



+ 



[v + l){2v + J + I] 
4 



E 



Ja+M-l 



•E 



> + l)(z/ + 2) 



E 



+ 



— (Et 

2c •''"2 



2Et _ , +Et _ ) + 



1 /2i/ + J - 1 



C2 



E^ 



,v~1 



{2u + J)E+^_ 



2z/+ J + 1 



E 



Jcr- 



Oil 
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and 



V X E+ _ 

Ja V 



E 



1 + J)(z/ + J + 1) 



Jo- 1/ ' ^2 



E 



+ j+l)(2^ + j+1) 

+ ^Ja-,v-\^ 



6z/2 + 6 Jz/ + 12z/ + J2 + 5 J + 6 



E 



Ja^ ,v — 2 ' ^ 

(i^ + l)(2z/ + J + 3) 



4 



(.H-l)(.H-2) 



./cr- 



-[(z/ + J)(j. + J + 1)E+^^^ + 2(zy + l)(z/ + J + 1)E+^^^^^, + + l)(^ + 2)E%,^^^,] 



1 f{u + J){u + J+l){2u+r 



2c2 

(z/+ l)(z/ + 2)(2z/ + 2J + 3)^, 

-E 



■E 



+ J+l)(2z/ + J + 2)E+^+^ 



+1 



JcT+,U + 2 



o 



For Ej^^, the right hand side of Eq. (l49l) and Eq. fl50|) must be multiphed by —1 and 
Ej^^ must be replaced everywhere by Ej^^. 



8. Outlook 



The results obtained here are an excellent starting point for calculating corrections to the 
paraxial approximation. Due to the equivalence to lowest order in ^ of the Laguerre- 
Gaussian solutions to the paraxial approximation and the spheroidal wave functions 
considered here, solutions to physical problems based on the paraxial approximation 
are effectively lowest-order solutions in | using spheroidal coordinates. Higher-order 
terms in | can then be included as perturbations to obtain results to any desired degree 
of accuracy, provided c is large enough that the resulting series converge. This procedure 
can be used in particular to calculate corrections to the eigenfrequencies and eigenmodes 
of a Fabry-Perot resonator. The resulting corrections can be observed experimentally 
in a high- finesse resonator [5]. 

A similar procedure as the one used here to write the radial functions in a new 
form could possibly be applied for solutions to the wave equation in other coordinate 
systems, in particular for elliptic coordinates. Such solutions should reduce to the 
Hermite-Gaussian solutions to the paraxial wave equation in the short wavelength limit. 
Applying such solutions to calculate eigenfrequencies and eigenmodes of two-dimensional 
resonators would allow comparison with results obtained previously for two-dimensional 
resonators using other methods [151 [HI flT] . 
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Appendix A. Recursion relations for Laguerre polynomials [13] 

L^-\x) = L^r'\x) - LtV\x) (A.l) 
xL^J^\x) = -(z/ + m)4"!i(^) + (2z/ + m + " + 1)4+1 (^) (A.2) 

= (i. + m)4™-i)(a;) - (i^ + l)^"^'^^;) (A.3) 

^ A^M (a;) = (a;) _ + m)LS (x) (A.4) 
ax 

Appendix B. r,n^ for p = 3 and p = A 

(3) - J_ /^ 4(;/ + !)(;/ + 2)(;/ + 3)(t^ + m + l)(;/ + m + 2)(;/ + m + 3) 
""""^^^ " 192 + 
4:/(i^ - - 2)(i/ + m)(j/ + Tn-l)(jy + m- 2) 

6(i/+ + 2)(i/ + m + + m + 2){iy'^ + mv - Av - 2m - 2) 

Qv{v - l){v + m)(:/ + m - + mi/ + + 3m + 3) 

3(z/ + + m + l)(14z/3 + 21mi/2 - 13i/2 + Tm^jy - 13mi^ - 2Av - Am^ - 12m - . 
Zv{v + m)(14i^3 + 21mj/2 + 55j/^ + 7m^i/ + 55mi/ + 44i/ + llm^ + 22m + 11) \ 



(B.l) 



(4) . , _ J_ / 2(;y + l)(t/ + 2){v + 3)(;/ + 4){v + m + 1)(;/ + m + 2)(;/ + m + 3)(;/ + m + 4) 

2v{v - l)(i/ - 2)(i^ - 'i){v + m)(j/ + m - l)(i^ + m - 2)(i^ + m - 3) 

4(i/ + l){v + 2){v + 3)(i^ + m + l)(i^ + m + 2)(i^ + m + 3)(i/^ + mi^ - 6;^ - 3m - 3) 

41^(2^ - l){v -2){v + m){v + m - l)(i/ + m - 2)(t/2 + ynzy + + 4?n + 4) 

(i^ + 1)(;^ + 2){v + m + 1)(;/ + m + 2) 
(l + ze)2 

{v^ + 2mv^ - 54i^3 + m^i^^ _ 8i„i,^2 ^ ^^^^2 _ 27^2^ + lllmi/ + I8O1/ + 30m^ + 90m + 60)+ (B.2) 
vlv — + m){i' + m — 1) 

{1^^ + 2mv^ + 58;^^ + mV^ + 87mz^^ + 279;^^ + 2%rr?v + 279mi/ + 208i/ + 58m,^ + 104m + 46)+ 



+ '^)^+^^^ + (-j,6 _^ (-3^^ _^ 2)i/5 + (3^2 + 5m - 311)1/"* + (m^ + 4m2 - 622m + M)v^+ 
(m^ - 382m2 + 126m + 653)1^^ + (-71m^ + 102m2 + 653m + 540)i/ + 30m-'^ + 162m2 + 270m + 13 

^^^^^^ (j,6 ^ (3^^ ^ 4)^5 ^ /3^^2 _^ -^Q^^ _ 306);/'' + (m^ + 8m2 - 612m - \?,2%)v^ + 2m^z/2- 
1 — z4 ^ 

(376?«2 + 1992m + 1470);/^ - (70m,^ + 866m2 + 1470m + 734)^ - (lOlm'^ + 323m2 + 367m + 145)) 
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